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ABSTRACT 

Competitive accretion, a process to explain the origin of the IMF, occurs when stars in 
a common gravitational potential accrete from a distributed gaseous component. Stars 
located near the centre of the potential benefit from the gravitational attraction of the 
full potential and accrete at much higher rates than do isolated stars. We show that 
concerns recently raised on the efficiency of competitive accretion are incorrect as they 
use globally averaged properties which are inappropriate for the detailed physics of a 
forming stellar cluster. A full treatment requires a realistic treatment of the cluster 
potential, the distribution of turbulent velocities and gas densities. Accreting gas does 
not travel at the global virial velocity of the system due to the velocity-sizescale 
relation inherent in turbulent gas and due to the lower velocity dispersion of small-N 
clusters in which much of the accretion occurs. Accretion occurs due to the effect of 
the local potential in tunneling gas down to the centre. Stars located in the gas-rich 
centres of such systems initially accrete from low relative velocity gas attaining larger 
masses before needing to accrete the higher velocity gas. Stars not in the centres of 
such potentials, or that enter the cluster later when the velocity dispersion is higher, 
do not accrete significantly and thus retain their low-masses. In competitive accretion, 
most stars do not continue to accrete significantly such that their masses are set from 
the fragmentation process. It is the few stars which continue to accrete that become 
higher-mass stars. Competitive accretion is therefore likely to be responsible for the 
formation of higher-mass stars and can explain the mass distribution, mass segregation 
and binary frequency of these stars. Global kinematics of competitive accretion models 
include large-scale mass infall, with mean inflow velocities of order « 0.5 km/s at 
scales of 0.5 pc, but infall signatures are likely to be confused by the large tangential 
velocities and the velocity dispersion present. Finally, we discuss potential limitations 
of competitive accretion and conclude that competitive accretion is currently the most 
likely model for the origin of the high-mass end of the IMF. 

Key words: stars: formation ~ stars: luminosity function, mass function - globular 
clusters and associations: general. 



1 INTRODUCTION 

One of the primary goals of a complete theory for star for- 
mation is to explain the origin of the distribution of stellar 
masses, the initial mass function (IMF). There have been 
many theories for the IMF that have involved a variety of 
physical processes from fragmentation (Zinnecker 1984; Lar- 
son 1985, Elmegreen 1997, Klessen, Burkert & Bate 1998; 
Klessen 2001;Bate, Bonnell & Bromm 2003; Bonnell, Clarke 
& Bate 2006), turbulence (Elmegreen 1993; Padoan & Nord- 
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lund 2002) accretion (Zinnecker 1982, Larson 1992, Bonnell 
et al. 1997,2QQlb; Klessen & Burkert 2000; Bate & Bon- 
nell 2005), feedback (Silk 1995, Adams & Fatuzzo 1996), 
magnetic fields (Shu. Li & Allen 2004) and a combination of 
these (Adams & Fatuzzo 1996). The simplest possible mech- 
anism that can explain the origin of the IMF is one that 
relies primarily on the physics of gravity, namely fragmen- 
tation, accretion and dynamical interactions. (Bonnell, Lar- 
son & Zinnecker 2006). In this scenario, gravitational frag- 
mentation of molecular clouds sets the mean stellar mass at 
~ O.5M0. Lower-mass stars and brown dwarfs arise due to 
the small Jeans masses produced in collapsing regions (fila- 
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ments, discs), followed by stellar interactions and potentially 
ejections to stop any subsequent accretion (Bate, Bonnell & 
Bromm 2002a; Reipurth & Clarke 2001). Higher-mass stars 
form due to continued accretion in a clustered environment 
where the overall system potential funnels gas down to the 
centre of the potential, to be accreted by the proto-massive 
stars located there. Not only can this reproduce the stellar 
IMF, but it also is able to account for the mass segregation 
of young stellar clusters and the binary properties of low 
and high-mass stars (Bonnell & Bate 2005; Bate, Bonnell & 
Bromm 2002b). 

Turbulent fragmentation has also been suggested as 
an alternative mechanism whereby the masses of stars arc 
determined by the size and densities of turbulent shocks 
(Padoan & Nordlund 2002). Lower velocity shocks produce 
weak but large shocks and hence result in higher mass 
clumps/stars whereas high-velocity shocks produce strong 
but narrow shocks and hence lower mass clumps/stars. One 
difflculty with this model (see also Ballesteros-Paredes et 
al.2006) is that the massive stars should form well sepa- 
rated and hence in relative isolation, not in the centre of 
dense stellar clusters. Indeed, in turbulent fragmentation, it 
is the lower-mass cores that are expected to form in such 
densely packed regions. 

Competitive accretion relies on the inefficiency of frag- 
mentation such that there is a large common reservoir of 
gas from which the protostars can accrete (e.g. Bonnell 
et al. 2001a). Observations of pre-stellar structures and of 
young stellar clusters both support this view with the large 
majority of the total mass being in a distributed gaseous 
form (Motte et al. 1998; Johnstone et al. 2000, 2004; Lada 
& Lada 2003; Bastian & Goodwin 2006). The second re- 
quirement is that the gas be free to move under the same 
gravitational acceleration as the stars. If the gas is fixed 
in place due to magnetic fields then accretion will be lim- 
ited. When these two requirements are filled, the dynam- 
ical timescale for accretion and evolution are similar such 
that a significant amount of gas can be accreted. Models of 
the formation of a stellar cluster show that the initial frag- 
mentation produces objects of order the Jeans mass of the 
cloud, even in the presence of turbulent velocities (Bonnell, 
Vine & Bate 2004). Accretion of nearby low relative veloc- 
ity material within their tidal radii increase these masses 
somewhat and helps produce the shallow IMF for low-mass 
stars (Bonnell et al.2001b; Klessen & Burkert 2000). The 
stars then fall together to form small stellar clusters which 
grow by accreting stars and gas (Bonnell et al.2003, 2004). 
This produces a strong correlation between the number of 
stars in the cluster and the mass of the most massive star 
it contains, as is found in young stellar clusters (Weidner 
& Kroupa 2006). Stars near the centre of the clusters bene- 
fit from the full potential which funnels gas down to them, 
increasing the local gas density. Once the stars establish a 
stellar dominated region in the centre of the clusters, they 
virialise and have higher relative gas velocities such that the 
Bondi-Hoyle radius, based on local parameters, now deter- 
mines the accretion rate. Accretion in this regime produces 
the higher-mass stars (Bonnell et al.2004) as well as the 
Salpeter-hke high-mass IMF (Bonnell et al.2001b; Bonnell 
et al.2003). Furthermore, the stellar dynamics maintain the 
accreting massive star in the centre of the system, and hence 



with a low relative velocity compared to the cluster poten- 
tial. 

Recently, Krumholz, Klein & McKee (2005a) have cast 
some doubt on this process by claiming that, accretion in 
such environments cannot significantly increase a star's mass 
and therefore does not play a role in establishing the stellar 
IMF. In part this is correct as with competitive accretion, 
most stars do not continue to accrete. It is the few that do 
that are important in terms of forming higher mass stars 
and the IMF. We have reanalysed our results in view of 
their work in order to show why accretion does occur to 
form higher-mass stars in our simulations and to establish 
if this is a realistic outcome of star formaiton. The principle 
difference is that the Krumholz et al. analysis used global 
parameters for the gas properties which can be significantly 
different from the local properties of a forming stellar clus- 
ter. We find that using global properties significantly under- 
estimates the correct accretion rates for the few stars that 
achieve higher masses. In §2 we discuss the Krumholz et al. 
analysis and point out some limitations in their approach. In 
§3 we reanalyse our numerical simulations. In §4 we discuss 
observational constraints of the models and in §5 we investi- 
gate potential limitations of competitive accretion. Finally, 
we summarize our current understanding of competitive ac- 
cretion in S6. 



2 ANALYTICAL ESTIMATES OF THE 
ACCRETION RATES 

The Krumholz et al. (2005a) approach considers accretion 
in a molecular cloud supported against self-gravity by its 
internal turbulent motions. In order to make the problem 
analytically tractable, they use globally averaged properties 
of the gas and relative velocities. In order to do this, they 
first assume that the gas is not globally self-gravitating. This 
implies that that the gas mass is then less than the mass 
in stars, and under such circumstances it is obvious that 
accretion cannot significantly alter a stars mass. They also 
neglect the larger scale cluster potential as being anything 
but a boundary to keep the gas within the system. Together, 
these assumptions ignore the similar accelerations that both 
the gas and star undergo in a cluster potential, and which 
to a large degree determine the relative velocities. 

The accretion rates can be estimated from the Bondi- 
Hoyle formalism where the accretion rate, M* is given by 



M* 



47rp 



{GM,) 



(1) 



where M* is the stellar mass, p the gas density, G the gravi- 
tational constant and v the relative velocity of the gas. Using 
such a formalism thus depends on having accurate values for 
the gas density, stellar masses and relative velocities. Unfor- 
tunately, Krumholz et al. (2005a) use globally averaged esti- 
mates for the gas density and the velocity instead of the local 
variables as used in the Bonnell et al. (2001a, 2004) numer- 
ical simulations. This causes the many order of magnitude 
difference in the calculated accretion rates. 

A stellar cluster or other self-gravitating system is not 
of uniform density but has a significantly increased density 
near the centre of the gravitational potential. For example, 
the mass density in the centre of the Orion nebula Clus- 
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Table 1. Accretion rates for global and local properties 



Property 


Global 


Local 


Gas density (g cm~"^) 


2x 10^19 


10-17 


Velocity dispersion (km s~^) 


2 


0.4 


Stellar mass (Mq) 


0.1 


0.5 


Accretion rate (Mq yr"-') 


10-9 


10-4 



ter (ONC) is approximately 50 times the mean mass den- 
sity of the system (Hillenbrand & Hartmann 1998). Sim- 
ilar central condensations are produced in simulations of 
turbulent fragmentation and cluster formation (Bonnell et 
al. 2003). Indeed, in competing models of massive star for- 
mation through the pressure- induced compression of turbu- 
lent cores, the central mass density is much higher than this 
(McKee & Tan 2003). Using a global virial velocity for the 
relative gas velocity is also problematic as it neglects that 
turbulence has a velocity size-scale relation and thus much 
local gas travels at significantly lower velocities. Further- 
more, if a core is formed through turbulent shocks, then 
much of the kinetic energy in the surrounding gas will also 
have been dissipated. The velocities of the forming stars are 
also not large for the same reason and even once they fall 
into local potential minima to form small-N clusters, the ve- 
locity dispersion in these systems is low. A forming cluster 
of 10 stars in a 0.1 pc radius will have a velocity dispersion 
of order 0.4 km/s instead of 2 km/s for the ONC. In terms 
of initial masses, Krumholz et al. use M* — O.IM0 rather 
than in the simulations where the mean initial masses are 
M^, ~ 0.5Mq. Thus, in terms of a Bondi-Hoyle accretion 
yields 



M, 



Af, 



Pglob V Af,,glob 



^'glob 



(2) 



M.,giob 

which when we incorporate the differences between local and 
global values gives 



= 50 X 5^ X 5^ 



M, 



1.5 X W 



(3) 



glob 



This difference of 1 x 10^ in the accretion rates is more than 
sufficient to explain why the Krumholz et al. analytical anal- 
ysis significantly underestimated the accretion rates from 
the numerical simulations. Thus, for example, using a rel- 
ative gas velocity of « = 0.4 km/s, a stellar mass of 
M, = O.5M0 and a gas density of p = 10" ^'' g cm"^ (^ 
the mass density in the core of the ONC), the accretion rate 
is then M» « 1O~*M0 yr""*^; whereas with v = 2 km/s, 
M* = O.IM0 and p = 2 x lO"^'' g cm"^, in line with 
Krumholz et al. (2005a), the estimated accretion rate is then 
M* « 1O~^M0 yr"!. These values are summarised in Ta- 
ble 1. Such a large difference in accretion rates for seemingly 
small changes in the physical parameters highlights the dan- 
ger of using global properties to predict what is, after all, a 
very local physical process. It also highlights the power of 
competitive accretion as small changes in the local proper- 
ties such as an acceleration from a stellar encounter and even 
the ejection into a lower gas-density region, can effectively 
halt the continued accretion. 

We can estimate the timescale for accretion to form a 
massive star by considering that a star forms with masses 
typical of the local Jeans mass which is near the median 




Mass 



Figure 1. The accretion timescale in years, from equation (8), 
is plotted as a function of the final stellar mass in M0. The plot 
assumes an initial mass of O.5M0 and a gas density of p = 10" ^^^ 
g cm— 3. The accretion timescale is directly proportional to the 
inverse of the stellar density. 



stellar mass (Af, ~ O.5Af0). We further consider that the 
relative gas velocity is initially low as expressed above (0.4 
km/s) and that it follows a turbulent velocity sizescale rela- 
tion of the form (Larson 1981; Heyer & Brunt 2004) 

v(xR^^^. (4) 

Assuming the mass distribution follows that of a centrally 
condensed p oc R~^ sphere such that M oc R (see Fig. 4 of 
Bonnell et al.2004), provides a relation between the turbu- 
lent velocity and the mass of the region considered 

V oc Af^/^ 



Combining this with equation (1) above gives 

dM, , {GM,f 



dt 



■' 4-Kp- 



(«i)' 



- M, \3/2' 



(5) 



(6) 



which can be integrated to yield an accretion timescale to 
form a star of final mass Af, of 



2vt m: 



tm — 



,1/2 



Af. 



1/2 



4npM,i G^ 



3/2 



(7) 



For typical values found in the numerical simulations of ui « 
0.5 km/s, pciuB « 10"^'' g cm~^ and M,i « O.5Af0, this 
reduces to 



tM ~ 10* 



0.5Mq 



1/2 



-1 



years. 



(8) 



This accretion timescale is plotted as a function of the final 
stellar mass in Figure 1. Thus to form a WMq star requires 
~ 3.3 X lO'* years while a 5OAf0 star requires ~ 9 x 10'' 
years. These are certainly reasonable numbers suggesting 
that competitive accretion can indeed account for the forma- 
tion of higher mass stars. There is still the issue of radiation 
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Figure 2. The 3-D velocity dispersion of the accreted gas is plotted (left) as a function of the cumulative accreted gas mass (in Mq) 
for four massive stars. The velocity dispersion is low for low masses due to the nature of turbulence and the low velocity dispersions in 
small-N clusters where much of the accretion occurs. The velocity dispersion increases with accreted mass due to the velocity-sizoscale 
relation for turbulent clouds. The accretion timescale, iacc ~ Mt/M*, is plotted as a function of accreted mass (right) showing that 
accretion occurs on relatively short timescales and thus determines the final masses of high-mass stars. 



pressure impeding the formation of stars greater than 10-40 
M0(Yorke & Kriigel 1979; Wolfire & Casinelh 1986), but 
sufficient current hypotheses exist to overcome this obsta- 
cle, including disc accretion and radiation beaming (Yorke 
& Sonnhalter 2002), Rayleigh- Taylor instabilities in the ac- 
cretion flow (Krumholz et al. 2005b) and stellar collisions 
(Bonnell, Bate & Zinnecker 1998; Bonnell & Bate 2002) and 
binary mergers (Bonnell & Bate 2005). 



3 NUMERICAL MODELS OF COMPETITIVE 
ACCRETION 

One of the advantages of using the smoothed particle hydro- 
dynamic (SPH) technique to follow star formation is that it 
is Lagrangian and therefore allows for the direct tracing of 
the fluid flow. We can thus directly probe the physical con- 
ditions of the accretion and assess what is occurring. In this 
section, we re-analyse the results from our previous simu- 
lation of the fragmentation of a turbulent cloud and the 
formation of a stellar cluster (Bonnell et al. 2003). It is 
this simulation which showed how massive stars can form 
through competitive accretion in the centre of the forming 
clusters (BonneU et al. 2004). 

The SPH technique uses sink-particles (Bate et al. 1995) 
to follow the formation of stars and any continuing gas ac- 
cretion onto them by adding the mass (and linear/angular 
momentum) of accreted particles to that of the sink-particle. 
It is therefore straightforward to track the evolution of the 
gas before it was accreted. In this way, we can determine 
the physical properties of the gas that formed a given star 
at all stages before it is accreted. For example, the physical 
distribution of this mass is shown in Bonnell et al. (2004) . 
Here, we analyse the velocity dispersion of the gas in order 
to determine how the accretion proceeds. 



In Figure 2 we plot the 3-D velocity dispersion, at the 
time the sink-particle forms, of the gas which eventually 
comprises four massive stars. This shows the velocity dis- 
persion of all the gas which is yet to accrete onto the star 
at the time the star first forms, and is plotted as a function 
of the mass in order of its accretion history. This therefore 
neglects any subsequent decrease in the gas's kinetic energy 
due to shocks. We see from figure 2 that the velocity disper- 
sion is flat at « 0.5 km/s for low-masses until M* « O.5M0. 
This is approximately the mean fragment mass (O.5M0) in 
the turbulent simulation and the velocity dispersion can be 
understood as being due to the gravitational collapse of the 
fragment at the sizescale of the fragment of several hundred 
AU. Accretion and higher-mass star formation accompanies 
the formation of a stellar cluster (Bonnell et al. 2004), such 
that the growth from low-masses starts in small-N clusters 
where the velocity dispersion is intrinsically low. Beyond 
this mass, the velocity dispersion increases with mass, as 
one expects for a turbulent medium, until nearly reaching 
the system velocity dispersion of 2 km/s. The actual in- 
crease in the velocity dispersion approximately follows the 
V tx M^", as discussed above. This is in keeping with a 
turbulent medium where the gas kinematics have not been 
overly affected by any dissipation, and further reenforces our 
choice of the initial velocity dispersion above of 0.5 km/s. 

We can use the above velocity dispersion as a function of 
accreted gas mass to evaluate typical accretion timescales. 
Figure 2 also shows the accretion timescale for a star to 
accrete its own mass from 



ta 



M, 

'mI'' 



(9) 



where M« is calculated from equation (1) using the veloc- 
ity dispersion and accreted mass from the simulations and 
assuming a gas density of lO"^'^ g cm— 3. At any given 
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time during the simulation after star formation has com- 
menced, at least 50 — IOOMq of gas has densities greater 
than this value. We can see from Figure 2 that the accre- 
tion timescale is between 10* to 10^ years for masses greater 
than the initial fragment mass of ~ 0.5Mq. From equa- 
tion (8), the estimated timescale for a star to double its 
mass is face ~ 1-4 x 10* years, on par with some of the 
smaller values for face in Figure 2. These timescale are suffi- 
ciently short that there is no difficulty forming higher mass 
stars from continued accretion in a clustered environment. 
The simulation formed two « 30Mq stars and a number of 
m > IOM0 stars in a few xlO^ years. The accretion rates 
(oc M^ /v^) actually increase with stellar mass, even though 
the gas now comes with larger velocity dispersions. This is 
due to the v oc M^ " seen in Figure 2. We can therefore con- 
clude that it is essential to use local quantities to correctly 
evaluate the accretion rates. 




4 INFALL AND LARGE-SCALE KINEMATICS 

Krumholz et al. (2005a) have argued that one way that com- 
petitive accretion could occur is if the system was signifi- 
cantly subvirial such that the relative gas velocities are then 
very small, and the whole system is in a state of radial col- 
lapse. We have seen above that this is not necessary due to 
the velocity-sizescale relation observed in turbulent molec- 
ular clouds such that smaller scales naturally have lower 
velocity dispersions. Furthermore, as the initial burst of ac- 
cretion occurs in small-N clusters where the stellar velocity 
dispersion is low, accreting stars can grow to higher-masses 
before needing to accrete higher velocity dispersion gas. Nev- 
ertheless, there must be some infiow into the central regions 
in order for the mass to arrive at where the massive star is 
forming. This, in fact is not unique to competitive accretion 
models as any model of massive star formation in the centre 
of a stellar cluster must at some point have inflow of mass 
to the cluster centre. In the McKee & Tan (2003) model for 
example, this inflow is simply presumed to have occurred at 
an earlier stage. 

We have analysed the gas kinematics from a simulation 
of a forming stellar cluster (Bonnell et al.2003). We find, in 
contradiction to the prediction of Krumholz et al. (2005a), 
that the system maintains approximate global virial equi- 
librium (ie, Ovir = i5kin/|Spot| ~ 0.5). This can be easily 
understood in the terms of gravitational dynamics where 
even a cold coUisionless system forms a virialised system at 
half its initial radius. In terms of a forming stellar cluster, as 
long as the gas is highly structured, and especially as it con- 
tains significant initial tangential motions, then it too will 
maintain a near-virial configuration. The kinetic energy of 
the system is initially equal to the potential energy. The sub- 
sequent decay of kinetic energy in shocks allows the system 
to become globally bound. Once gravity forms local poten- 
tial wells into which the stars and gas fall, the gravitational 
acceleration increases both the radial and the tangential ve- 
locities such that the gas has larger, and still highly disor- 
dered motions. In order to attain the evolution suggested 
in Krumholz et al.(2005a), there would have to be no initial 
tangential motions or no conservation of angular momentum 
during infall. 

In order to quantify the kinematics, we plot in Figure 3 



Figure 3. The mean radial velocity (km/s) in sliells is plotted 
as a function of the distance (in pc) from the forming massive 
star (symbols). The dashed lines denote the velocity dispersion 
around the mean radial velocity while the solid line is the mean 
tangential velocity. Note that the inward radial velocity is smaller 
than both the tangential velocity and the dispersion in the radial 
velocity. 



the mean radial and tangential velocities as a function of 
radius from the accreting massive star. The SPH particles 
are binned as a function of radial shells and the mean ra- 
dial and tangential velocities are calculated along with the 
dispersion in radial velocity. The mean radial velocity is in- 
wards from ~ 0.6 pc with a typical value of 0.5 — 1 km/s. 
The velocities outside this radius are generally outwards due 
to the large kinetic energy in the initial conditions. The dis- 
persion in the radial velocity is actually larger (±1 km/s) 
than this inward velocity indicating that the motions are 
very chaotic. The tangential velocities are also larger than 
the radial velocity demonstrating that the system is not un- 
dergoing a simple collapse process but involves significant 
chaotic/turbulent motions and is kinematically hot. The sys- 
tem is significantly different from the cold collapse model 
suggested by Krumholz et al. (2005a). 

Observational studies have searched for infall signa- 
tures in massive star formation. Recently, Fuller et al.(2005), 
Motte et al. (2005) and Perreto et al.(2006) have detected of 
such motions through asymmetric line profiles, and by spa- 
tially resolved kinematics, respectively. The estimated infall 
velocities are of the order 0.1 to 2 km/s, in rough agree- 
ment with the mean inward velocities seen in the numerical 
simulations. 

Observed linewidths in regions of massive star forma- 
tion have been found to be of order a few km/s (e.g. Garay 
2005; Beuther et al.2006), This is significantly larger than 
the values discussed above in terms of the local velocity dis- 
persion of the gas on relatively small scales. The difference 
can be understood as being due to the large lengthscale of 
the core probed along the line of sight. Even when small 
areas of the core are observed, the gas is likely to lie along 
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Figure 4. The mass of the gas within a projected radius of 0.1 
pc from the forming massive star is plotted as a function of its 
line-of-sight velocity at two different times during the evolution 
when the most massive star has 4 and 12 Mq, respectively. The 
mass is plotted in units of MQpc~^(km/s)~-'^ while the velocities 
are in units of km/s. Although local velocity dispersions are small 
(si 0.5 km/s Figure 2), the line of sight velocity dispersions are 
significantly larger ( few km/s ) due to the multiple clumps in the 
line of sight. 



an extended. This effect shown in Figure 4 where we plot 
the hne of sight velocity distribution of the gas mass within 
a projected areas of radius 0.1 pc centred on the forming 
massive star. This effective velocity profile for an optically 
thin tracer has a line width of several km/s in agreement 
with observations. Thus, although the small-scale velocity 
dispersion in a 3-D volume can be fairly low, the line of 
sight velocity dispersion in a 2-D projected area is much 
larger due to the greater sizescale probed. The existence of 
multiple distinct clumps can also be seen in the right hand 
panel where two clumps are now kinematically distinct. Such 
highly structured profiles which can be observed along some 
lines of sight are due to the hierarchical fragmentation of 
the turbulent cloud (Bonnell ct al. 2003). 



5 POTENTIAL LIMITATIONS OF 
COMPETITIVE ACCRETION 

We have seen above that competitive accretion arises natu- 
rally from the fragmentation of a turbulent molecular cloud 
where there is a common (cluster) potential such that stars 
near the centre of the potential accrete from higher density 
gas and generally have lower relative velocities. The relative 
gas velocities are low due to the nature of turbulence and 
its velocity-sizescale relation with lower velocities on smaller 
sizescales. In order for competitive accretion not to work as 
suggested by Krumholz et al. (2005a), requires that the local 
gas velocities are very high, generally as high as the turbu- 
lent velocity on the largest scale of the cloud. Neglecting that 
this would violate the velocity-sizescale relation {v ex R^'^), 
we discuss here potential mechanisms that could decouple 
the gas and stellar kinematics. 

In order to completely decouple the stellar and gas kine- 
matics requires a driving source for the gas motions. Such 
a driving source, potentially related to maintaining the tur- 
bulent support for many dynamical lifetimes, is often spec- 
ulated to be the internal feedback from low-mass stars. The 
problem is that although the energetics from low-mass stel- 
lar jets and outflows are sufficient to offset the decay of tur- 



bulent energy, it is very unclear if they can input sufficient 
kinetic energy into the system without completely remov- 
ing the gas (Arce et al.2006). Jets are commonly found to 
escape the bounds of their natal clouds suggesting that the 
majority of their momenta is deposited at large distances, 
outside the main star forming region Stanke et al.2000; see 
also Arce et al.2006; Bally et al.2006). Even the intrinsically 
isotropic feedback from high-mass stars escapes in prefer- 
ential directions due to the non-uniform gas distributions 
(Dale et al. 2005; Krumholz et al. 2005c). The feedback 
decreases the accretion rates but does not halt accretion. 
Indeed, in the toy-model feedback simulation by Li & Naka- 
mura (2005), the energy injection is sufficient to disrupt the 
core but is quickly damped from the system allowing it to 
continue its unimpeded collapse. 

A further problem is one of balancing the timescales. 
Energy injection from jets and outflows is relatively quick 
with typical dynamical times of the flows is generally of order 
10^ to lO"* years (e.g. Reipurth & Bally 2001; Beuther et al. 
2002). Giant outflows have longer dynamical timescales but 
considering outflow velocities of 10 to 100 km/s, the energy 
injection timescale of an outflow before it leaves he limited 
size of a forming stellar cluster is at most a few W^years. 
This timescale is smaller than the dynamical timescale of the 
system as a whole (few 10^ years). It is difficult to envision 
how the feedback can be exactly tuned to support the system 
without disrupting it. If excessive energy is deposited by the 
feedback, there is no opportunity for the cluster to adapt 
before the gas is dispersed. Likewise, if insufficient energy is 
deposited, then the system will continue to form stars. As 
the timescale for star formation is similar to the dynamical 
timescale of the cluster, it cannot be halted on the dynamical 
timescales of the outflows such that feedback is highly likely 
to overshoot the energy balance and completely unbind the 
system. Indeed, observations suggest that this is common 
(Arce et al.2006) and that there is actually no need for long 
lifetimes as star formation appears to generally occur on 
dynamical timescales (Elmegreen 2000). 



6 CONCLUSIONS 

We have shown that concerns raised about the efficacy of 
competitive accretion are misplaced due to the misuse of 
global variables to estimate what is a local physical process. 
In competitive accretion, it is the few stars located early on 
in the centre of local stellar systems that become higher- 
mass stars. Competitive accretion starts with the accretion 
of low relative velocity gas due to the nature of turbulence 
and the low velocity dispersion in small-N clusters. Accre- 
tion proceeds to higher relative velocity gas once the star has 
attained sufficient mass to maintain a high accretion rate. 
Stars that enter a stellar cluster later or that reside far from 
the centre of the potential do not accrete significantly and 
hence lose the competition to accrete from the distributed 
gas. In this way, competitive accretion sets the distribution 
of stellar masses by determining which stars accrete to at- 
tain higher-masses and which do not so that they remain 
lower-mass stars. Although infall does occur, as it needs to 
in any model for massive star formation in the centre of 
a cluster, emission line signatures are likely to be confused 
due to the large chaotic and tangential motions present. The 
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forming cluster never approaches the status of a cold, radi- 
ally collapsing system. 



ACKNOWLEDGMENTS 

We thank Henrik Beuther, Cathie Clarke, Ralf Klessen, 
Richard Larson, Jim Pringle and Hans Zinnecker for dis- 
cussions and comments which helped formulate this letter. 
We acknowledge the contribution of the U.K. Astrophysi- 
cal Fluids Facility (UKAFF) in providing a facility to test 
and refine our models. MRB is grateful for the support of a 
Philip Leverhulme Prize. 



REFERENCES 

Adams F. C, Fatuzzo M., 1996, ApJ, 464, 256 
Arce H. G., Sheperd D., Gueth F., Lee C-F., Bachiller R., Rosen 
A., Beuther H., 2006, in Protostars and Planets V, eds B. 
Reipurth et al., in press 
Ballesteros-Paredes, J., Klessen R., Mac Low M.-M., Vazqucz- 
Semadeni E., 2006, in Protostars and Planets V, eds B. 
Reipurth et al., in press 
Bally, J., Reipurth B., Davis C. J., 2006, in Protostars and Planets 

V, eds B. Reipurth et al., in press 
Bastian N., Goodwin S. P., 2006, MNRAS, in press. 
Bate M. R., Bonncll I. A., 2005, MNRAS, 356,1201 

Bromm V., 2002a, MNRAS, 332, 65 
Bromm V., 2002b, MNRAS, 336, 705 
Bromm V., 2003, MNRAS, 339, 577 
Price N. M., 1995, MNRAS, 277, 362. 
B., mcKee C. F.,Tan J. C., 2006, in 
Protostars and Planets V, eds B. Reipurth et al., in press 
Beuther H., Schiike P., Sridharan T. K., Mcnten K. M., Walmsley 
0. M., Wyrowski P., 2002, A&A, 383, 892 



Bate M. R., 


Bonnell I. A., 


Bate M. R., 


Bonnell I. A., 


Bate M. R., 


Bonnell I. A., 


Bate M. R., 


Bonnell L A., 


Beuther H., 


Churchwell E, 



Bonnell I. A., 
Bonnell I. A., 
Bonncll I. A., 

RAS, 285 
Bonnell I. A. 

MNRAS, 
Bonnell I. A., 
Bonnell LA., 

93 
Bonnell I. A., 
Bonncll I. A. 

MNRAS, 
Bonnell LA., 



Bate M. 

Bate M. 

Bate M. 
, 201 

Bate M. 
323, 785 

Bate M. 

Bate M. 



R., 2002, MNRAS, 336, 659 
R., 2005, MNRAS, 362, 915 
R., Clarke C. J., Pringle J. E., 



1997, MN- 



R., Clarke C. J., & Pringle J. E., 2001a, 



Vine S. G., 2003, 
& Zinnecker, H., 



MNRAS, 343, 413 
1998, MNRAS, 298, 



Clarke C. J., Bate M. R., 2006, MNRAS, in press. 
, Clarke C. J., Bate M. R., Pringle J. E., 2001b, 
324, 573. 

Larson R.B., Zinnecker H., 2006, in Protostars and 
Planets V, eds B. Reipurth et al., in press 
Bonnell I. A., Vine S. G., Bate M. R., 2004, MNRAS, 349, 735 
Dale J. E., Bonnell I. A., Clarke C. J., Bate M. R., 2005, MNRAS, 
358, 291. 

ApJ, 419, L29. 

ApJ, 486, 944. 

ApJ, 530, 277 

MNRAS, 354, 367 

S. J., Sridharan T. K., 2005, A&A, 442, 



G., 


1993, 


G., 


1997, 


G., 


2000, 


G., 


2004, 


Wi 


Uiams 



Elmegreen B. 

Elmegreen B. 

Elmegreen B. 

Elmegreen B. 

Fuller G. A., 
949 

Garay G., 2005, in Massive star birth, a crossroads of astro- 
physics, lAU Symp 227, eds R. Cesaroni et al., CUP, p. 86. 

Heyer M. H., Brunt C. M., 2004, ApJL, 615, L45. 

Hillenbrand L. A., Hartmann L. W., 1998, ApJ, 492, 540 

Johnstone D., Wilson C. D., Moriarty-Schieven G., Joncas G., 
Smith G., Grcgerscn E., Fich M., 2000, ApJ, 545, 327. 

Johnstone D., Di Francesco J., Kirk H., 2004, ApJ, 611, L45. 

Klessen R. S., 2001, ApJ, 556, 837 



Klessen R. S., Burkert A., 2000, ApJS, 128, 287 

Klessen R. S., Burkert A., Bate M. R., 1998, ApJ, 501, L205 

Krumholz M. R., McKee C. F., Klein R. L, 2005a, Nature, 438, 

332. 
Krumholz M. R., Klein R. L, McKee, C. F., 2005b, in Massive 
star birth: A crossroads of astrophysics: lAU Symposium 227, 
eds R. Cesaroni et al., (CUP), 231. 
Krumholz M. R., McKee C. F., Klein R. I., 2005c, ApJL,618, L33. 
Lada C. J., Lada, E. 2003, ARA&A, 41, 57 
Larson R. B., 1981, MNRAS, 194, 809. 

214, 379 
256, 641. 
359, 211 

astro, arXiv:astro-ph/0512278 
ApJ, 585, 850 



Larson R. B., 1985, MNRAS, 

Larson R. B., 1992, MNRAS, 

Larson R. B., 2005, MNRAS, 

Li Z.-Y., Nakamura P., 2005, 

McKee, C.F., Tan, J.C, 2003, 

Motte P., Andre P., Neri R.,1998, A&A, 336, 150. 

Motte F., Bontcmps S., Schiike P., Lis D. C, Schneider N., 

Mcnten K. M., 2005, lAUS, 227, 151 
Padoan P., Nordlund A., 2002, ApJ, 576, 870. 
Padoan P., Nordlund A., 2004) ApJ, 617, 559. 
Peretto N., Andre P., Belloche A., 2006, A&A, 445, 979 
Reipurth B., Clarke C, 2001, AJ, 122, 432. 
Reipurth B., Bally, J., 2001, ARA&A, 39, 403 
Shu F. H., Li Z.-Y., Allen A., 2004, ApJ, 601, 930. 
Silk J., 1995, ApJ, 438, L41 

eidner C, Kroupa P., 2006, MNRAS, 365, 1333 
Wolfire M.G., CassineUi J. P., 1987, ApJ, 319, 850 
Yorke H., Kriigel E., 1977, A&A, 54, 183 
Yorke H., Sonnhalter, C, 2002, ApJ, 569, 846 
Zinnecker H., 1982 in Symposium on the Orion Nebula to Honour 

Henry Draper, eds A. E. Glassgold et al.. New York Academy 

of Sciences, 395, p. 226 
Zinnecker H., 1984, MNRAS, 210, 43 



